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ABSTRACT 


-  I 

This  report  is  intended  as  a  modern  sequel  to  'Band  Pass  Filters  with  Linear 
Phased,  published  over  20  years  ago.  Here,  the  linear  phase  algorithms  are  adapted  to 
low-pass  integrated  active  circuit  topologies  for  potential  use  as  the  input  filter  in 
sampled  data  control  systems  and  other  low-pass  applications  for  which  traditional 
coil-capacitor-transformer  filter  designs  are  no  longer  appropriate.  _ 

The  linear  phase  algorithm  is  applicable  to  a  variety  of  circuit  topologies  which 
are  not  mathematically  equivalent.  Application  to  the  overall  transfer  function  (the 
"simple?  approach)  results  in  transmission  that  is  nearly  equal-ripple  in  both  amplitude 
and  phase,  the  errors  being  equivalent  to  echos  in  a  misterminated  transmission  line. 
The  resulting  filters  compare  favorably  with  modest  order  Tchebychef  filters  (lacking 
phase  control).  Application  to  branches  of  a  feedback  topology  (mathematically 
equivalent  to  a  lattice)  leads  to  control  of  echoes  and  stop-band  zeros  at  the  cost  of 
extra  elements.  In  this  case,  the  results  compare  favorably  with  ordinary  Butterworth 
filters  (lacking  phase  control)  having  the  same  total  number  of  poles.  „ 

Practical  means  for  inserting  stop  band  zeros  without  redesign  of  the  pass-band  are 
considered.  Circuits  of  both  the  R-C  and  switchcap  types  are  discussed 
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Section  1 

INTRODUCTION 


Some  years  ago,  the  writer  described  algorithms  for  designing  slab-sided  flat- 
topped  band-pass  Alters  whose  in-band  phase  response  is  an  inherently  linear  function 
of  frequencyfl].  Elsewhere,  he  also  disclosed  adaptations  of  the  algorithms  to  band¬ 
pass  Alters  with  quadratic  phase[2]  ("chirp"  Alters)  and  to  non-integer  power  law 
impedances[3].  Colleagues  independently  adapted  the  linear  phase  technique  to  sam¬ 
pled  data  for  use  in  Vocoders[4]  and  general  digital  Alter  design[$]. 

An  application  to  the  design  of  low-pass  Alters  was  promised  [1]  but  was  not  pub¬ 
lished.  A  vigorous  sampled  data  technology  suggests  examination  of  the  low-pass  lin¬ 
ear  phase  Alter  as  a  means  of  preAltering  such  data  before,  or  concurrently  with,  sam¬ 
pling.  Today,  a  high-Q  single  pole-pair  is  easily  synthesized  with  op-amps,  resistors, 
and  capacitors,  either  as  an  analog  feed-back  network[6]  [7],  or  as  a  switched  capaci¬ 
tor  network[8].  Such  responses  are  easily  joined  with  isolation  ampliAers  and  summing 
nodes,  all  on  a  single  integrated  chip.  This  report  emphasizes  topologies  appropriate  to 
such  active  integrated  circuit  synthesis. 

Figure  1  shows  a  typical  amplitude  template  for  specifying  a  single  stage  low  pass 
linear  phase  Alter.  Here,  M  is  the  number  of  "main  poles".  A  linear  phase  design  also 
requires  a  minimum  overhead  of  two  "corrector  poles". 

The  template  differs  from  ordinary  speciAcation  templates  in  (a)  allowing  a  wider 
transition  between  the  1  dB  point  and  the  beginning  of  stop  band  roll-off,  and  (b) 
allowing  an  additional  transition  between  the  tolerance  limits  deep  in  band  and  the 
1  dB  point.  The  template  includes  the  possibility  a  stop  band  zero  for  supressing  car¬ 
rier  leakage.  In-band,  the  phase  is  bounded  by  the  same  template,  with  0.1  radian 
phase  distortion  replacing  the  1  dB  amplitude  bound. 

Linear  Phase  (L-Ph)  designs  compare  favorably  efficiency  in  use  of  numbers  of 
elements  with  Butterworth  Alters  for  given  tolerances  of  pass-band  transmission  and 
rapidity  of  reaching  Axed  stop  band  attenuation.  Figure  2  shows  a  nearly  maximally 
Aat  L-Ph  Alter  transmission  having  13  total  poles  in  the  s-plane  (11  so-called  main 
poles,  a  pair  of  "correctors,"  and  a  pair  of  "compensators").  Its  transmission  compares 
favorably  with  that  of  a  Butterworth  having  17  total  poles.  The  linearity  of  the  phase, 
within  a  few  degrees  of  a  pure  delay  over  most  of  the  pass-band  is  thus  "for  free". 
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Figure  1.  Template  applicable  to  single  stage  linear  phase  filter. 
Stop  band  slope  of  template  is  12.SM  dB  per  u0. 


Some  of  the  results  reported  below  appear  to  be  new.  These  include  the  single 
and  dual  loop  realizations  of  L-Ph  lattice  networks  (Sections  2  and  8),  details  of 
manipulation  of  band-edge  behavior  (Section  S)  and  out-of-band  zeros  (Section  7). 
In  addition,  many  results  of  computer  simulation  are  presented,  along  with  design 
equations  which  the  reader  can  use  to  meet  his  own  particular  set  of  requirements. 

There  remains  a  widespread  impression  that  phase  linearity  and  pass-band  square¬ 
ness  cannot  be  designed  together.  This  impression  is  the  legacy  of  an  era  in  which  fil¬ 
ters  were  used  almost  exclusively  for  frequency  division  telephony,  for  which  the  noto¬ 
rious  phase  insensitivity  of  the  ear  led  designers  to  concentrate  on  amplitude  only. 
Even  so,  in  a  1931  patent[9]  Bode  came  very  close  to  finding  the  L-Ph  algorithms 
used  here  (he  used  pole-zero  instead  of  pole-residue  synthesis).  The  present  writer 
believes  that  if  serious  commercial  applications  had  then  existed,  Bode  would  have 
succeeded.  It  remains  to  be  seen  whether  there  are  sufficient  commercial  applications 
to  result  in  the  widespread  application  of  these  techniques  today. 
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Figure  2.  Comparison  of  Linear  Phase  with  Butterworth  designs. 
See  Figure  13  for  further  details. 
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Section  2 


NETWORK  TOPOLOGY 


Assume  a  mathematical  description  of  a  network  as  a  rational  transfer  function 
H(w),  often  specified  by  its  zeros  and  poles.  The  L-Ph  design  algorithms  all  specify 
the  poles  pk  and  the  residues  rk  of  the  equivalent  partial  fraction  expansion  of  H(w) 
[1][8]: 


H(w)  =  c  +  r,/(w-p,)  +  r2/(w-p2)  +  ...  +  rm/(w-pm)  (1) 

The  L-Ph  algorithm  was  introduced  in  the  s-plane.  Gold  and  Rader  [5],  and  recently 
Lish  [8],  have  recognized  that  if  the  algorithm  is  applied  in  the  z-plane,  a  sampled 
data  L-Ph  Alter  also  results  (although  it  is  not  exactly  the  same  filter).  The  w  in 
Eq  (1)  is  intended  to  represent  either  the  Laplace  complex  frequency  s,  a+ju,  or  the 
sampled  data  transform  variable  z. 


Figure  3  shows  an  active  circuit  version  of  the  L-Ph  Alter  topology  based  on  the 
use  of  H(w)  as  a  Alter  response.  1  This  will  be  called  the  ’  simple"  L-Ph  Alter  topology. 
Each  block  represents  an  active  circuit  having  the  labelled  transmission.  Signal  flow  is 
in  the  direction  of  the  arrows. 


The  simple  Alters  were  called  "inAnite  loss"  Alters  in  the  original  paper  [1]. 


Physical  realizability  of  H(w)  requires  that  the  individual  terms  occur  in  pairs 
having  complex  conjugate  pole  positions  and  residues.  Each  such  pair  is  first  summed 
to  produce  a  physical  2nd  order  system  yr  (A  pole  on  the  real  w-axis  remains  a  1st 
order  system.)  Those  y,  of  H(w)  whose  real  part  is  positive  near  d.c.  are  driven  from 
the  "+"  output  of  the  input  isolation  stage.  The  sum  will  be  called  "YA".  Those  terms 
whose  real  part  is  negative  near  d.c.  are  driven  from  the  side;  their  sum  will  be 
called  "Yb".  The  individual  transmissions  are  simply  added  to  produce  the  output. 

H(w)  =  Ya(w)-Yb(w)  (2) 

When  the  poles  and  residues  of  H(w)  are  chosen  in  accordance  with  the  L-Ph  algo¬ 
rithm  below,  the  transmission  of  such  a  filter  behaves  like  a  finite-bandwidth  section 
of  transmission  line  with  internal  reflections  from  mistermination  at  the  ends.  These 
echos  result  in  correlated  nearly-equal  ripple  errors  in  both  amplitude  and  phase 
within  the  pass  band. 
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Figure  4.  Realizations  of  the  H  (w)  topology. 

(a)  "Low  Loss"  topology,  (b)  "Feedback"  topology. 


as 


In  addition  to  "simple"  filters,  the  original  paper  explored  the  use  of  Y  A  and  YB  « 
branch  admittances  feeding  a  load  admittance  G,  as  shown  in  Figure  4a.  The  corre¬ 
sponding  active  network  will  be  called  here  the  "feedback"  or  "loop"  topology.  When 
this  new  network  is  driven  by  a  voltage  source,  the  transmission  is  given  by  H  (w) 


H(w)  =  (Ya-Yb)/(Ya+Yb+G) 


O) 


Of  course,  this  is  not  H(w).  However,  if  Y  .  and  YB  have  been  chosen  according  to 
the  L-Ph  algorithm,  the  resulting  filter  is  L-Ph.  Moreover,  there  is  is  a  choice  of  G 
that  terminates  the  equivalent  delay  cable  and  thus  supresses  the  echos. 


In  general,  G  is  the  sum  of  a  constant  G0  and  a  frequency  dependent  G  p(w).  It  is 
evident  from  Eq  (3)  that  H(w)  will  not  change  if  only  G0  is  assigned  to  the  amplifier 
while  Gp/2  is  added  to  each  of  YA  and  YB.  Such  an  assignment  is  convenient  when 
G  is  also  given  in  pole-residue  form. 

The  original  paper  also  considered  a  third  generic  network,  having  finite  source 
impedance.  Such  a  network  is  shown  in  circuit  form  on  the  left  hand  side  of  Figure  S. 
One  reduction  of  this  case  to  the  feedback  topology  is  shown  in  Figure  3.  Two  grand 
feedback  loops  are  now  required,  and  two  extra  frequency-dependent  gain  elements  G. 
The  Figure  includes  a  reduction  of  G  to  the  form  of  a  fixed  gain  G0  and  feedback  ele¬ 
ments  G  . 

p 

The  transmission  of  this  new  topology  is  given  by  H(w) 


Figure  3.  Matched  source  topology. 


H(w)  =  2(  Y  .-Y  _)/ (G+ 2  Y  A)(G+ 2  Y  _) 


Both  one  loop  and  two  loop  topologies  are  treated  further  in  Sections  6,  7,  and  8 
below. 
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Section  3 


LINEAR  PHASE  ALGORITHM 
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(a)  ODD  POLE  PATTERN 


(b)  EVEN  POLE  PATTERN 


Figure  6.  Allowable  patterns  for  low  pass  poles. 


The  L-Ph  design  algorithm  places  all  of  the  "main"  poles  at  equal  intervals  2 a 
along  a  line  parallel  to  the  ju-axis  in  the  s-plane,  located  at  a  distance  b  behind  that 
axis.  The  residues  at  these  poles  all  have  the  same  magnitude;  but  residues  at  adjacent 
poles  have  alternating  signs.  Beyond  upper  and  lower  band  edges,  the  line  of  main 
poles  is  simply  stopped. 

Nominal  band  edges  at  u=±<jo  are  taken  to  lie  half  way  between  the  last  main  and 
the  first  out-of-band  pole  of  the  periodic  array.  If  there  are  M  main  poles  (counting 
both  half  planes),  then 

o0  =  Ma  or  a  =  u«/M  (5) 


Let  E(ju)  be  the  difference  between  the  transmission  of  the  main  poles  alone,  and 
that  of  an  infinite  chain  of  such  poles  multiplied  by  an  ideal  slab-sided  amplitude.  In 
[1]  it  was  shown  that  E(ju)  can  be  approximated  by  yet  another  pole  located  at  each 
band  edge,  with  a  residue  equal  to  one-half  that  of  the  first  adjacent  out-of-band  pole 
dropped.  These  are  the  corrector  poles.  The  nominal  locations  are  on  the  same  line 
as  the  main  poles,  at  s=±jw0+b. 


For  band-pass  filters,  the  entire  s-plane  design  is  carried  out  in  the  upper  half¬ 
plane.  Physical  reality  requires  another  set  of  poles  in  the  lower  half  plane.  This 
merely  adds  another  pass  band  at  negative  frequencies,  whose  skirts  are  too  low  to 
affect  the  transmission  of  the  design  filter  (and  vice  versa).  At  low  pass,  the  physical 
requirement  of  conjugate  residues  at  conjugate  pole  positions  can  conflict  with  the 
L-Ph  requirement  for  alternating  signs  at  adjacent  pole  positions.  Only  two  patterns  of 
low-pass  pole  positions  satisfy  both  sets  of  requirements:  an  even  pole  number  pattern 
and  an  odd  number  pattern.  Both  are  shown  in  Figure  6. 

The  usual  pattern  is  that  with  an  odd  number  of  poles,  which  includes  one  pole  on 
the  real  axis.  All  of  the  examples  in  this  report  are  calculated  on  the  assumption  that 
the  number  of  poles  is  odd.  The  even-number  pattern  is  useful  mainly  as  a  means  of 
inducing  ripple  reduction  in  a  cascade  with  a  simple  filter  having  an  odd-number  pole 
pattern.  The  even  pattern  is  not  applicable  to  loop  topologies. 

In  all  these  formulations,  there  is  an  arbitrary  scale  factor,  which  has  been  set 
equal  to  b  below.  A  single  contribution  to  H(s)  is  thus  y 

ym  =  (-j)”  b/(s-jma+b)  (6) 

where  m  is  an  even  integer  between  ±M  for  odd  pole  filters.  The  corrector  residue  is 


rc  =  0.5  (~j)M+l  b 


(7) 


Figure  7.  Position  of  low  pass  poles  in  z-plane. 


Figure  7  shows  the  location  of  poles  of  a  simple  L-Ph  filter  specified  directly  in 
the  z-plane.  In  this  case,  the  main  poles  of  H(z)  are  y  m(z)  They  are  located  at  equal 
intervals  2a  along  a  circle  whose  center  is  at  z  = 0 . 
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yn(z)  =  (-j)“  />/(z-e,“#+/>)  (8) 

where  p<  1  is  the  radius  of  the  pole  circle.  Lish  1  noted  the  explicit  transformation 
between  such  a  a  z-plane  design  and  an  s-plane  design: 

z  =  e*r  or  s  =  (1/r)  log(z)  (9) 

where  r  is  the  sampling  interval  in  the  time  domain. 

However,  the  simple  s-plane  and  z-plane  designs  are  not  the  same  filter.  The 
s-plane  design  that  would  correspond  to  Lish’s  simple  z-plane  filter  has  periodic  pass 
bands  at  intervals  of  1/r  in  radian  frequency.  This  periodicity  of  pass  bands  will  be 
dealt  with  later,  in  Section  10. 


Section  4 

SIMPLE  FILTER 


The  L-Ph  algorithm  defines  a  pass  band  H(w)  that  is  nearly  periodic. 


H(s)  = 

I(-l)k  b/(s+b-2jka) 

for  odd  pole  filters 

H(s)  = 

I(-l)k  jb/(s+b-j(2k+  l)a) 

for  even  pole  filters 

(10) 

If  there  are  no 

limits  on  the  k’s,  Eqs  (10)  can  be  summed  to  a  closed  form[10]. 

H(s)  = 

irb/[2a  sinh(n(s+b)/2a)] 

[odd] 

(a) 

H(s)  = 

trb/[2a  cosh(n(s+b)/2a)] 

[even] 

(b) 

(11) 


From  these  equations  it  follows  that  the  maxima  H  ^  and  minima  Hmta  of 
transmission  are 

H„„  =  trb/2a  sinh(nb/2a) 

Hmjn  =  nb/2a  cosh(nb/2a)  (12) 

It  also  follows  that  the  ratio  of  the  squared  magnitudes  of  maximum  to  minimum 
periodic  transmission  is  1  +€*,  where 

€  =  l/sinh(rrb/2a) 

•X.  2  e-’b/2*  (13) 

The  last  two  equations  state  that  low  ripple  is  accompanied  by  low  absolute  transmis¬ 
sion  resulting  from  the  small  differences  between  overlapping  y  ..  Thus  the  design  is 
subject  to  a  tradeoff  between  low  theoretical  ripple  and  sensitivity  to  errors  in  the 
individual  yr  To  aid  the  selection  of  b/a,  Figure  8  gives  the  average  magnitude  |H(j  u)| 
and  the  amplitude  of  the  fundamental  component  of  the  ripple  as  functions  of  b/a. 

Figure  9  exhibits  typical  simple  L-Ph  filter  responses  with  nominal  correctors. 
That  transmission  does  not  quite  fit  the  template  in  Figure  1.  This  can  be  cured  by 
minor  adjustments  to  the  corrector,  discussed  in  the  next  section.  The  Figure  includes 
the  result  of  one  such  adjustment;  another  is  shown  in  Figure  11.  Such  corrector 
modifications  have  limited  scope.  For  applications  which  require  robust  suppression 
of  apparent  echos,  carrier  leakage,  or  out  of  band  interference,  it  is  necessary  to  use 
cascades  of  filters  and/or  filters  in  the  loop  configurations. 


Figure  9.  Response  of  filter  near  band  edge. 

Nominal  corrector  position  and  residue. 

Corrector  modified  to  95%  residue,  95%  b,  and  position  0.1a  out-of-band  also  shown. 


4.1  COMPARISON  WITH  ORDINARY  TCHEBYCHEF 

• 

It  is  natural  to  compare  the  transmission  of  the  nearly  equi-ripple  simple  L-Ph 
filter  to  that  of  the  standard  equi-ripple  Tchebychef  for  phase  distortion  in-band  and 
attenuation  slope  at  the  onset  of  the  stop  band.  The  comparison  is  made  difficult  by 
differences  in  the  scaling  rules  for  the  design  of  the  two  filter  types. 

Except  when  M  is  small  enough  to  produce  special  cases  (M<5),  the  phase  distor¬ 
tion  of  the  L-Ph  design  occurs  within  1.5u0/M  of  the  band  edge  at  u0,  while  the  stop 
band  behavior  is  generally  a  function  of  >? 


rj  =  M(u-u0)Ajo 


(14) 


In  particular,  the  nominal  template  stop  band  attenuation,  Figure  1,  is  given,  near  the 
band  edge,2  by  W (>?) 


Observe  that  this  is  not  the  result  of  analysis  of  a  given  filter,  but  rather  a  specifica¬ 
tion  slope. 
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On  the  other  hand,  the  transmission  T(j  u)  of  an  N-pole  Tchebychef  is  given  by 

|T(ju)|l  =  (1  +  e1  cosh1  Ncosh-'(l+r>/M)l-'  (16) 

Evidently,  this  transmission  depends  on  e,  M,  and  N,  as  well  as  i?,  where  that  of  the 
L-Ph  can  be  specified  in  terms  of  n  only. 

To  cast  the  preceding  result  into  a  form  in  which  it  can  be  compared  directly  with 


Eq  (15),  one  can  make  use  of  the  approximation 3 

cosh”'(l+u)  =  (2u)w  +  . higher  order  terms  mu'4  (17) 

cosh'O+rj/M)  =  (In/ M)“  +  .  (18) 

to  obtain 

rr(jo)|3  -  (4A1)  exp[-2N (2r> /M)  Vi  ]  (19) 


Obviously,  this  result  and  Eq  (15)  are  not  comparable.  For  any  £  it  is  always  possible 
to  find  an  N  such  that  the  Tchebychef  attenuation  lies  either  above  or  below  the 
M-pole  L-Ph  for  the  first  50  dB  of  attenuation,  and  vice  versa. 

Figure  9  also  shows  the  stop  band  transmission  of  two  standard  Tchebychefs  hav¬ 
ing  the  same  in-band  ripple  and  a  slightly  narrower  nominal  bandwidth.  Curves  are 
plotted  on  the  assumption  that  M=N,  so  that  number  of  L-Ph  main  poles  equals  the 
number  of  Tchebychef  poles.  The  L-Ph  stop  band  attenuation  shown  is  bounded  by 
that  of  5-  and  8-pole  Tchebychefs  out  to  a  point  at  which  the  L-Ph  filter  "bottoms 
out"  (either  by  design  or  on  account  of  random  construction  errors),  other  situations. 
If  the  b/a  in  the  figure  had  been  1.25,  the  bounding  Tchebychefs  occur  at  8-  and 
11-pole;  if  b/a  had  been  1.50,  the  bounding  Tchebychefs  occur  at  11-  and  15-pole. 

The  phase  distortion  of  a  Tchebychef  is  simpler  than  its  attenuation.  It  oscillates 
around  an  overall  distortion  curve  A 

A  =  M  [arcsin(u/u0)  -  u/u0J  (20) 

The  onset  of  unacceptable  phase  distortion  thus  occurs  deeper  into  the  pass-band  as 
N  is  increased.  If  for  comparison  with  L-Ph,  the  edge  of  unacceptable  phase  distortion 
is  taken  to  be  at  0.1  radian  distortion,  then  this  edge  occurs  at  u/u0  of  roughly  0.5  for 
a  5-pole  and  roughly  0.4  for  a  8-pole  Tchebychef.  In  contrast,  in  the  L-Ph  filter,  as 
the  number  of  main  poles  M  is  increased  the  span  of  acceptable  phase  distortion 
moves  closer  to  the  band  edge;  for  a  0.1  radian  limit,  the  phase  distortion  is  satisfac¬ 
tory  over  (M-l)/M  of  w«. 


Higher  order  terms  in  the  expansion  can  be  obtained  using  the  identity 

coshM(l+y)  =  sinh',((2y+y,}V4) 

together  with  the  power  series  expansion 

sinh''(x)  =  x  -  x5/6  +  3xs/40  - . 


Section  S 


MODIFICATIONS  OF  SIMPLE  FILTERS 


This  section  considers  several  methods  of  tailoring  the  response  of  simple  filters  by 
adjustment  of  the  corrector  poles  and  by  use  of  higher  order  correction. 

The  correctors  y  c  are  an  approximation  to  the  ideal  difference  E(j  o).  Figure  10 
shows  the  real  and  imaginary  parts  of  typical  E(j  u),  and  of  corresponding  nominal 
correctors  Y^.4  The  left  hand  side  of  the  Figure  illustrates  b/a  less  than  one.  The 
right  hand  side  is  typical  of  b/a>l.  (The  right  hand  amplitude  scales  have  been 
expanded  to  show  the  approximation  error  on  the  same  relative  scale  as  for  b/a=0.5). 
In  both  cases,  the  missing  halves  of  the  curves  are  obtained  by  reflection  with  respect 
to  the  band  edge,  with  the  imaginary  parts  changing  sign.  In  the  1930s  Bode[ll]  gave 
relationships  that  hold  along  s=ju  between  the  real  and  imaginary  parts  of  physically 
realizable  H(s).  Key  concepts  included  minimum  phase  for  a  given  gain  curve,  maxi¬ 
mum  loss  for  a  given  phase  curve,  minimum  reactance  for  a  given  real  part,  and  vari¬ 
ous  Hilbert  Transform  5  relationships  between  real  and  imaginary  parts  when  one  or 
the  other  is  specified  for  all  frequency. 

Although  Re(E)  and  Im(E)  have  the  correct  general  relationship  for  physical  real¬ 
izability,  E(ju)  is  not  "minimum  reactance".  It  follows  that  there  is  no  physical  "super- 
corrector"  arbitrarily  close  to  E(ju)  for  all  frequencies.  The  curve  labelled  "HT"  shows 
the  minimum  reactance  imaginary  part  for  Re(E)  when  b/a  =0.3,  obtained  from  Re(E) 
by  numerical  Hilbert  transform. 

However,  excellent  correction  can  be  obtained  outside  a  transition  zone  straddling 
the  band  edge.  In  this  report,  the  width  of  that  zone  is  arbitrarily  fixed  to  be  ±a 
around  u0. 

Adjustments  to  the  corrector  pole(s)  of  a  simple  filter  can  effect  modest  tradeoffs 
between  extending  L-Ph  into  the  transition  zone  (at  the  expense  of  increased  attenua¬ 
tion  in-band)  versus  extending  the  flat-topped  amplitude  response  towards  the  band 
edge  (at  the  cost  of  increased  phase  distortion).  Most  of  the  illustrative  examples 
involve  such  adjustments  to  single  pole  correctors. 


4  A  filter  with  an  odd  number  of  poles  is  assumed.  If  the  number  of  poles  is  actually 
even,  statements  about  "real"  and  "imaginary"  parts  must  be  interchanged.  Such  an 
interchange  has  no  effect  on  the  overall  results. 

5  If  F(x)  is  a  suitably  bounded  generalized  function,  its  Hilbert  Transform  G(x)  is 
given  by 

G(x)  =  (1/rr)  /  F(u)du  /(x-u) 

where  the  integral  is  understood  in  the  Cauchy  (or  generalized  function)  sense  to 
run  between  u=-~  and  «. 


AMPLITUDE  —  UNITS  OF  b/a 


5.1  ADJUSTMENT  OF  SIMPLE  CORRECTOR  POLE 

The  following  observations  can  be  used  to  guide  adjustments  to  the  residue  or  position 
of  a  single  corrector  pole  to  make  minor  local  improvements  to  the  match  to  E(j  u). 

1.  Reducing  the  corrector  residue  from  nominal  (V£),  improves  the  match  of  both 
real  and  imaginary  parts  near  w0±a,  at  the  expense  of  a  broad  "bottom"  to  the 
in-band  ripple  and  stop  band  attenuation. 

2.  The  primary  effect  of  changing  corrector  b  is  to  change  the  real  part  of  yc. 
Thus  improving  the  knee  of  the  pass-band  at  the  zone  edge  generally  involves 
lowering  the  b  of  the  corrector  pole. 

3.  Moving  the  frequency  of  the  corrector  (parallel  to  the  j  u-axis  or  z-circle) 
improves  the  match  of  yc  and  E(u)  close  to  one  side  of  the  transition  zone,  at 
the  expense  of  the  match  on  the  other  side  of  the  zone. 

As  a  rule,  the  illustrative  examples  of  this  report  involve  such  adjustments  to  the  cor¬ 
rector.  For  example,  Figure  9  includes  modifications  to  the  corrector  pole  to  make  the 
in-band  response  conform  to  the  template  of  Figure  1. 


5.2  USE  OF  MULTIPLE  POLE  CORRECTORS 


Figure  1 1 .  Use  of  two  term  corrector. 


In  the  range  0.5<b/a<  1.5  two  poles  can  be  used  in  combination  to  produce  a  better 
match  (outside  the  zone)  to  E(ju)  than  is  available  from  a  single  pole.  Let  the  ampli¬ 
tudes  of  these  two  correctors  be  y ,  and  y2,  and  the  damping  b,  and  b2: 

yCi  =  i  A»*+b,) 

yC2  =  by2  /(s*+b2)  (21) 

with  the  signs  cf  the  y’s  chosen  to  give  the  same  sense  as  the  nominal  corrector.  The 
most  important  on  the  y’s  and  b’s  is  asymptotic  equality  of  Im(E)  and  Im(yc).  A  less 
important  constraint  is  the  asymptotic  equality  of  the  real  parts.  Together,  these  con¬ 
straints  suffice  to  determine  the  y’s  in  terms  of  the  b’s: 

y,  =  (b-b2)/(b,-b2) 

y2  =  (b— b , )/ (b2— b , )  (22) 

One  strategy  for  choosing  the  b’s  is  to  chose  one  corrector  pole  to  have  the  lowest  b, 
that  is  consistent  with  practical  hardware  and  boundedness  of  y  c  in  the  transition 


zone.  The  other  pole  is  taken  with  a  fairly  large  b .  (say  b/a-2.3).  The  gammas  are 
calculated  according  to  Eqs  (22).  Using  this  point  of  departure,  a  satisfactory  set  of 
b’s  and  y’s  is  found  by  computer  cut-and-try. 


TABLE  1 

Parameters  for  Two  Pole  Corrector 


b,/a 

bj/a 

>2 

Filter  b/a=0.5 

2.5 

0.09 

0.20 

0.41 

Eq(22) 

0.065 

0.435 

Filter  b/a  =1.0 

2.5 

0.17 

0.40 

0.33 

Eq(22) 

0.143 

0.357 

Filter  b/a  =1.5 

2.5 

0.25 

0.60 

0.25 

Eq(22) 

0.237 

0.263 

Figure  1 1  shows  the  use  of  two  pole  correctors  for  y  c  in  simple  L-Ph  filters  for 
b/a=  1  and  b/a=1.5.  The  actual  cut-and-try  parameters  used  are  given  in  Table  1. 
The  Table  also  lists  corresponding  y’s  computed  from  Eqs  (22).  Note  that  the  trans¬ 
mission  of  the  filter  for  b/a  =1  has  lower  distortion  in-band  and  sharper  attenuation 
skirts  than  attenuation  out-of-band  than  either  of  the  single  pole  adjustments  shown  in 
Figure  9. 

How  far  can  the  multiple  corrector  pole  approach  be  extended?  For  any  given 
topology  there  is  a  minimum  width  of  the  transition,  such  that  a  corrector  y  c  which 
matches  E(ju)  outside  this  zone  will  not  produce  unacceptable  peaks  in  transmission 
within  it.  Theoretical  specification  of  this  minimum  cannot  be  undertaken  here.  In 
fact,  for  low  b/a,  even  the  nominal  corrector  can  cause  unacceptable  peaks  or  "ears" 
inside  the  transition  zone,  especially  in  the  single  loop  feedback  topology  dicussed  in 
Section  6. 

Low  b/a  is  not  used  in  single  simple  filters.  However,  it  is  used  in  cascades  and  in 
loop  topologies.  Low  b/a  also  lays  bare  character  of  the  nominal  transition  zone 
edges.  For  b/a<<l,  the  |Im(E)|  tends  to  a  limit  of  .349b/a  at  ±a  from  the  band  edge, 
while  the  |Im(yc)  tends  towards  0.5b/a.  Both  real  parts  tend  to  vanish.  The  result  is 
a  net  residual  imaginary  part  ±0.1 5 lb/a  at  the  zone  edges.  This  can  be  compared  with 
a  periodic  pass  band  minimum  |H(u)|  of  nb/2a  as  b/a-0;  and  a  periodic  maximum 
Re(H)  of  unity. 

Thus,  the  residual  produces  transmission  at  roughly  -20  dB  with  respect  to  the 
pass  band  minima  at  the  outer  edge  of  the  transition  zone,  and  0. 15b/a  radian  phase 
distortion  at  the  inner  edge.  The  effect  of  these  residuals  can  be  reduced  in  general 
only  by  blunting  the  sharpness  of  the  transition  at  the  band  edge  by  the  use  of  high-b 
corrector  elements,  or  by  introducing  unacceptable  peaks  within  the  zone. 


Section  6 


FEEDBACK  TOPOLOGY 


IDEAL 

11  1  XFORMER 

Figure  12.  Lattice  network. 


Any  passive  bilateral  network  that  can  be  inserted  between  a  source  having  (com¬ 
plex)  admittance  G,  and  a  load  having  impedance  G2  can  be  realized  in  the  form  of 
a  half  lattice,  shown  in  Figure  12.  6  In  the  original  paper  it  was  shown  that  excellent 
L-Ph  filters  can  be  obtained  in  this  topology  by  applying  the  L-Ph  algorithm  to  Y  A 
and  Yb  rather  than  to  T(s).  This  use  of  the  algorithm  can  be  extended  to  low  pass 
Filters  with  an  odd  total  number  of  poles. 

Two  special  cases  are  of  interest:  the  case  of  "matched"  source  and  loads, 
G,=GI=G,  and  the  low  impedance  source,  G,=“,  G2=G.  For  matched  loads,  the 
transmission  has  already  been  given  as  Eq  (4).  In  the  low  impedance  case  (G  l=°°)  the 
current  source  is  replaced  by  a  voltage  source.  The  transmission  is  given  in  Eq  (3). 

T(s)  =  ( Y  A- Y  g)/[G+  Ya+ Yb]  (3) 


6  In  networks  textbooks,  main  attention  is  usually  given  to  is  a  full  symmetric  lattice 
with  balanced  inputs  and  outputs;  the  original  paper  made  use  of  full  lattice  analy¬ 
sis.  To  compare  equations,  the  half-lattice  G’s  in  this  report  are  twice  the  full  lat¬ 
tice  G’s  of  the  original  paper.  Branch  Y’s  are  unchanged  in  meaning. 


Let  yA  and  yB  be  the  infinite  pole-residue  arrays  that  are  eventually  truncated  into  the 
main  poles  of  and  YB.  The  echos  in  pass  band  transmission  are  eliminated  by 
chosing  G  according  to  the  image  parameter  rule 


G2  =  4yAyB  (23) 

In  particular,  if  the  total  number  of  main  poles  (after  truncation)  is  odd,  the  individ¬ 
ual  yA  and  yB  (as  well  as  y^-y^  are  periodic.  They  can  be  expressed  in  closed  form 
as  tangents.  For  example,  if  y  A  happens  to  include  the  pole  on  the  negative  real  axis: 

yA  =  ±(nb/4a)  tanh((s+b)/4a) 

yB  =  ±(rrb/4a)  coth((s+b)/4a)  (24) 

The  periodic  product  yAyB  is  independent  of  frequency.  The  image  admittance  G 
reduces  to  G0 

G„  =  nb/2a  (25) 

Direct  substitution  of  y  A,  yB,  and  G0  from  the  preceding  equations  into  Eqs  (3)  or  (4) 
produces  (after  some  complex  variable  trigonometry)  a  simple  delay  with  attenuation 
independent  of  frequency. 

T(s)  =  constant  x  e""b/1*  e~"'/2*  (26) 


6.1  SHUNT  COMPENSATION 

In  this  subsection  it  is  assumed  that  the  branch  YA  contains  the  main  pole(s) 
closest  to  the  band  edge  and  that  YB  contains  the  corrector  pole(s),  yc.7  Let  y*A  and 
y*B  represent  the  periodic  poles  that  are  dropped  when  the  network  is  truncated: 

Ya  =  yA  -  y*A  (main  poles)  (a) 

Yb  =  yB  -  y*B  +  yc  (main  poles  plus  corrector)  (b) 

G  =  G0  +  Gp  (c) 

(27) 

The  original  paper  points  out  that  the  choices 


7  If  Yb  contains  the  pole  closest  to  band  edge,  YA  includes  the  corrector,  and  Gp  is 
made  up  of  the  dropped  terms  y* 


yc  =  y*B  -y*A  (in  band  only) 


(a) 


=  (yA-yB>  -  (y*A-yV  (°ut  of  band) 

Gp  =  2y\  (b) 

(28) 

produce  ideal  pass-band  transmission  and  stop-band  attenuation.  The  reader  can  con¬ 
firm  this  by  using  Eq  (cpi)  to  reduce  Eq  (27),  then  substituting  the  result  into  either  of 
Eqs  (22).  The  actual  corrector  is  an  approximation  to  Eqs  (28a)  and  (28b).  Proper 
selection  of  the  corrector  is  all  that  is  needed  to  convert  the  filter  from  infinite  band¬ 
width  to  design  bandwidth. 

Eq  (27c)  introduces  a  second  kind  of  corrector,  called  in  the  original  paper  "shunt 
compensation,"  Gp.  Ideally,  this  new  Gp  is  the  sum  of  all  the  terms  y*  that  were 
dropped  from  yA  (but  not  yB)  when  it  was  truncated  from  periodic  to  finite  form.  7 

The  ideal  Gp  contains  the  poles  located  at  M+3a,  M+7a,  ....  and  at  -M-3a,  -M-7a . 

This  G  can  be  expressed  in  closed  form  using  the  (^-function  (the  logarithmic  deriva¬ 
tive  of  the  F  -function): 

Gp  =  (b/2aj)  foK CM+31/4  -(s+b)/4aj)  -  ^((M+3)/4  -Ks+b)/4aj)J  (29) 

In  practice,  shunt  compensation  consists  of  a  few  network  elements  which  approximate 
the  ideal  Gp  in-band.  Near  d.c.,  Gp  can  be  approximated  by  a  nearly  constant  real 
part,  and  an  imaginary  part  which  varies  linearly  with  frequency: 

Gp  =  go  b2/aJ  +  g,s  b/a  +  . 

g„  =  1/(M+  1)J  +  .... 

g,  =  1/(M+1)  -  (2/3+bV3aJ)/(M+l)]  +  .  (30) 

Near  the  band  edge,  the  simple  power  series  development  of  G  is  not  wholly  ade¬ 
quate.  In  addition,  the  synthesis  of  the  term  g  ,s  is  not  convenient  for  the  active  feed¬ 
back  topology.  It  is  better  to  use  an  out-of-band  pole  pair,  even  when  the  resulting 
zero  in  transmission  is  not  the  primary  objective. 

The  parameters  of  a  one  pole  pair  G  can  be  computed  using  Eq  (29).'  The  fre¬ 
quency  and  residue  are  chosen  so  that  In^Gp)  has  the  correct  value  at  the  in-band 
edge  of  the  transition  zone  (at  s*=ja),  and  the  correct  slope  at  d.c.  The  damping  of 


*  The  properties  of  the  {^-function  are  discussed  in  standard  treatments  of  higher 
transcendental  functions.  For  Re(w)>0,  <p(w)  is  ln(w)-l/2w-d?(w)/dw  where  f (w)  is 
a  continued  fraction.  (Handbook  of  Mathematical  Functions,  Abramowitz  and 
Stegun,  Editors,  Dover,  New  York,  1970,  §6.1.48).  In  reducing  formulas,  the  author 
has  used  the  recursion  (op  cit  §6.3.5): 
t//(z)  =  — 1/z  +  j^(z+l) 

an  asymptotic  expansion  for  dtft/dz  (from  op  cit  §6.4.12): 

f(z)  =  l/(z-V4)  -  l/24(z-l/2)5  +  . 

and  his  own  formula  for  numerical  computations: 
tp(x)  =  ln(x-.493)  +l/10x3  ±  0.0003  for  x  >  1.1 


the  pole  pair  is  next  selected  to  produce  the  proper  real  part  for  G  at  the  zone  edge. 
Finally,  g«  is  adjusted  so  that  the  real  part  is  correct  at  d.c. 

Let  p  and  q  be  the  normalized  magnitudes  of  the  real  and  imaginary  parts  of  EQ  u) 
at  the  zone  edge: 

E(ju0)  =  p  b'/aJ  -  jq  b/a  (31) 

Approximate  formulas  for  p  and  q  are* 
p  =  .198  +  1/4M  -  .026  bVa1 

q  =  -0.0579  +  ln(M)  -  .040  bVa  (32) 

The  result  of  the  matching  process  just  outlined  is  a  pole  pair  given  by  parameters 
(A,  B,  C): 

Locations  at  s  =  Bb  ±j  (M+A)a  Residues  =  Cb 

Table  2  gives  values  for  A,  B,  and  C  calculated  according  to  the  above  matching  pro¬ 
cess  for  the  case  b/a = Vi.  The  numbers  are  not  significantly  different  for  any  other 
b/a-0. 


TABLE  2 

Parameters  of  Nominal  Compensation.... b/a  ~Vi 


M 

q 

A 

A/q 

C 

B 

AG0 

5 

0.74 

6.52 

8.84 

10.74 

1.29 

0.24 

7 

0.91 

6.99 

7.72 

12.04 

1.23 

0.11 

9 

1.03 

7.64 

7.42 

13.72 

1.22 

0.07 

11 

1.13 

8.33 

7.36 

15.47 

1.23 

0.04 

13 

1.21 

9.01 

7.42 

17.23 

1.25 

0.03 

15 

1.29 

9.69 

7.53 

18.99 

1.28 

0.02 

17 

1.35 

10.35 

7.67 

20.73 

1.31 

0.02 

19 

1.40 

10.99 

7.83 

22.45 

1.34 

0.01 

21 

1.45 

11.63 

8.00 

24.17 

1.37 

0.01 

Figure  1 3  shows  the  transmission  of  a  single  loop  feedback  network  with  1 1  main 
poles  (5  each  in  upper  and  lower  half  planes)  using  compensation  from  Table  2  and 
nominal  correctors.  Figure  2  is  similar.  Here,  corresponding  curves  are  plotted  for  b/a 
equal  to  0.25  and  0.50.  Note  that  they  are  essentially  identical.  However,  at  b/a  =0.25 
there  is  a  serious  peak  or  "ear"  in  transmission,  which  gets  worse  as  b/a  is  reduced. 
This  peak  is  eliminated  in  the  "matched"  or  "double  loop"  topology. 


PHASE  ERROR  (deg)  MAGNITUDE  (dB) 


MAGNITUDE  (dB) 


|B(u)F  =  [1  +  (u/u.)‘N]" 

=  [1  +  exp  2Nln(u/u0)]-' 

=  [1  +  exp  2Nln(l+n/M)]- 

-  [1  +  exp  2N^/M]‘‘ 

-  exp(-2Nn/M)  (33) 

Evidently,  in  both  cases  the  attenuation  slope  follow  the  same  scaling  law.  Observe 
that  for  a  given  slope,  then  number  of  L-Ph  main  poles  is  substantially  less  than  the 
total  poles  of  the  Butterworth,  so  that  the  comparison  can  remain  favorable  to  L-Ph 
even  including  correctors  and  compensators. 

Figure  13  also  compares  the  transmission  and  phase  distortion  of  the  L-Ph  filter 
with  the  Butterworth  filter  whose  transmission  it  most  nearly  resembles.  This  Butter- 
worth  has  the  same  group  delay  at  d.c.  and  virtually  the  same  stop  band  transmission 
out  to  where  the  L-Ph  filter  "bottoms".  It  has  it  M/2  poles,  rounded  to  the  nearest 
integer.  Allowing  for  the  overhead  of  corrector  and  compensator,  there  are  still  just  as 
few  second  order  systems  in  the  L-Ph  filter! 

The  figure  also  shows  the  phase  distortion  of  the  Butterworth,  which  is  similar  to 
the  average  distortion  of  the  Tchebychef.  Note  that  it  exceeds  the  0. 1  radian  criterion 
near  midband.  Here,  the  comparison  favors  the  L-Ph  over  the  Butterworth  on  all 
practical  counts.  The  price  is  inexact  theoretical  flatness,  and  the  1  dB  shoulder 
within  2u0  of  the  band  edge. 
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Section  7 


MAGNITUDE  (dB) 


Figure  14  shows  filter  transmission  with  the  "ideal"  compensation  G  specified  by 
Eq  (28).  Since  YA  and  Y,  contain  no  poles  in  the  stop  band,  the  poles  of  G  (in 
accordance  with  Eq  (3))  insert  zeros  into  the  stop  band  transmission  at  regular  inter¬ 
vals  beyond  the  band  edge.  Compensation  that  was  concentrated  at  a  single  zero  in  the 
stop  band  is  now  spread  over  several  zeros.  Up  to  the  first  zero  the  attenuation 
exceeds  that  for  single  element  compensation.  Beyc  d  this  point  the  loss  proceeds  in 
stair-step  fashion  and  is  not  immediately  as  profound  as  for  the  single  element. 


Inasmuch  as  G0  is  also  a  part  of  G  in  the  stop  band,  it  follows  from  Eq  (3)  that 
the  attenuation  cannot  be  less  than  that  of  a  simple  filter  (Y  ^-Y^)/G0,  whatever  be 
the  form  of  G  ,  so  long  as  YA,  YB,  and  Gp  are  positive  real  (i.e.  individually  stable) 
elements.  This  minimum  attenuation  curve  is  also  shown  in  the  Figure.  Using  lower 
loss  elements  for  the  poles  of  Gp  will  result  in  deeper  stop  band  zeros;  but  the  max¬ 
ima  between  zeros  then  tend  to  touch  the  minimum  loss  curve.  There  is  also  an 
adjustment  to  G0  to  avoid  a  reduction  in  Re(G)  in  the  pass  band.) 


7.2  INSERTION  OF  ARBITRARY  STOP  BAND  ZEROS 

There  are  two  separate  stratagems  available  for  the  insertion  of  zeros  into  a  L-Ph 
filter  stop  band.  The  first  makes  use  of  adjustments  to  the  corrector;  it  is  available  in 
all  cases.  The  second  makes  use  of  the  compensation;  it  is  available  only  in  the  feed¬ 
back  topologies. 

It  is  always  possible  to  pick  the  damping  and  residue  of  the  corrector  in  such  a 
way  as  to  obtain  equality  of  E(j  u)  and  Yc  at  one  frequency  in  the  stop  band  (and  at  a 
corresponding  frequency  in  the  pass  band).  It  is  also  possible  to  achieve  the  same 
result  asymmetrically  by  adjusting  corrector  position  instead  of  corrector  damping. 
Such  zeros  are  available  in  all  L-Ph  topologies  because  they  affect  only  the  Y  A-Y„ 
portions  of  the  transmission  equations. 

However,  the  technique  is  sensitive  to  minor  changes  in  component  values.  It  is  not 
useful  for  maintaining  a  zero  against  fixed  carrier  frequency  leakage. 

The  second  stratagem  is  available  only  in  the  feedback  topologies.  It  is  to  make 
use  of  the  zeros  naturally  present  in  shunt  compensation.  The  procedure  has  three 
steps: 

1 .  First,  the  nominal  pole  of  G  closest  to  the  band  edge  is  deliberately  retained. 
The  remainder  of  Gp  is  approximated  by  a  single  second  order  element  (as 
before),  on  the  basis  of  the  remaining  G  '  at  the  zone  edge  and  the  remaining 
slope  g,'  at  d.c. 

2.  If  the  nominal  pole  pair,  call  it  G  J(  is  to  remain  in  its  natural  position  at 
s*=±j3a,  it  is  retained  with  its  full  residue  and  assigned  the  same  nominal  b  as 
a  main  pole.  However,  the  f>  may  be  lowered  to  achieve  a  deeper  zero. 

3.  The  frequency  of  G  3  can  be  moved  without  affecting  either  the  Re(G  J  or  the 
slope  of  ImfGp)  at  the  zone  edge  if  the  residue  and  the  damping  are  both 
adjusted  in  the  ratio  f 


f  =  ((<jj^0)/4a)J 


(34) 
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Section  8 

DOUBLE  FEEDBACK  CONFIGURATION 


The  transmission  shown  for  the  (single  loop)  feedback  configuration  of  Figure  13 
exhibits  an  "ear"  at  the  band  edge  for  b/a=0.25  in  the  single  loop  topology.  Such  ears 
appear  to  be  an  unavoidable  part  of  the  single  loop  transmission  for  b/a  less  than 
about  one-third. 

The  cure  utilized  in  the  original  paper  was  to  drive  the  network  from  a  finite 
impedance,  rather  than  zero  impedance  source.  This  approach  had  several  advantages: 

1.  A  low  loss  network  could  be  designed  to  operate  between  a  resistive  source 
and  a  resistive  load. 

2.  The  transmission  "ears"  go  away. 

3.  When  source  and  load  resistance  are  taken  equal,  the  pass-band  transmission 
is  insensitive  to  minor  variations  of  individual  components. 

4.  The  zeros  of  transmission  are  deeper  than  with  the  zero  impedance  source. 

Most  of  these  advantages  remain  when  the  mathematics  is  realized  in  the  two  loop 
topology  shown  in  Figure  3.  Figure  16,  for  example,  is  the  transmission  of  a  two  loop 
network  with  the  same  parameters  as  single  loop  Figure  13.  Note  that  the  peaks  in 
transmission  at  the  band  edge  are  gone. 

In  the  two  loop  topology  of  Figure  3,  the  elements  G  are  used  three  times.  As  a 
result 

1.  There  are  two  more  amplifiers  and  two  more  sets  of  elements  G  than  in  the 
single  loop  case. 

2.  An  unbalance  in  two  nominally  equal  branch  G  ’s,  actually  realized  along 
with  the  Y’s,  is  equivalent  to  an  adjustment  in  Yc.  This  can  be  an  extra 
degree  of  freedom  in  the  design  of  the  corrector,  or  source  of  error  in  realiza¬ 
tion. 

The  extra  depth  of  zeros  in  the  two  loop  topology  is  the  result  of  all  signal  paths 
having  to  pass  through  an  equivalent  "G"  element  twice.  This  advantage  is  accompa¬ 
nied  by  a  potential  6  dB  increase  in  the  transmission  at  other  stop  band  frequencies. 
This  increase  can  be  deduced  by  considering  the  ratio  9,,  of  the  two  loop  transmis¬ 
sion  (from  Eq  (4)  to  the  one  loop  transmission  (from  Eq  (3»,  normalized  so  that 
|©21|—  1  within  the  pass  band 


9  2t  =  2G0(YA+Y#)/(G+2YA)(G+2YB) 


(35) 
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Figure  16.  Two  loop  transmission. 


In  the  stop  band,  YA  and  YB  tend  to  zero,  while  away  from  its  poles  G  tends  to  G0. 
Thus,  asymptotically,  021  tends  to  2,  which  is  6  dB  in  power.  However,  much  of  this 
asymptotic  region  of  lower  attenuation  is  far  enough  away  from  the  band  edge  so  that 
theoretical  6  dB  has  little  practical  significance.  The  "matched  source"  configuration 
remains  the  topology  of  choice  if  lumped  passive  components  can  be  used.  9  The 
double  feedback  loop  also  appears  to  be  the  topology  of  choice  for  critical  active  cir¬ 
cuit  applications.  Otherwise,  the  single  loop  is  to  be  preferred  over  the  simple  topology 
for  most  active  circuit  realizations. 


9  The  main  pole  on  the  real  axis  contributes  an  inductor  to  the  branch  admittance, 
which  can  be  combined  with  the  ideal  transformer  to  produce  a  real  transformer. 
In  the  context  of  the  original  paper,  if  the  inductances  of  the  other  branches  y  ,  are 
all  L,  then  this  "d.c."  inductance  is  2L. 


Section  9 


TYPICAL  CIRCUITS 


Although  the  theoretical  network  building  block  is  the  individual  pole  and  residue,  the 
physical  building  block  is  the  sum  of  a  given  upper  half  plane  pole  and  its  complex 
conjugate  mate.  If  the  position  of  the  uhp  pole  is  at  s=jma+b,  then  the  physical 
building  block  is  yB(s)  if  the  number  of  main  poles  is  odd;  and  y  ^  if  M  is  even. 


y»o(s>  =  b(s+b)/(s2  +m2a2  +b2  +2bs) 

(•) 

y^s)  =  jmab/(s2  +m2a2  +b2  +2bs) 

(b) 

(36) 

Note  that  the  frequency  s=jma  is  not  a  good  place  to  align  either  of  these  2nd  order 
systems.  For  M  odd,  a  good  alignment  frequency  is 

umo  =  (mJaJ-bJ)“  (37) 

where  the  phase  and  magnitude  of  the  transmission  are  zero  degrees  and  unity,  respec¬ 
tively.  For  M  even,  a  good  alignment  frequency  is 

=  (m2a2  +b2)*  (38) 

where  the  phase  and  transmission  are  90  *  and  m/^(m2+b2),  respectively. 


Figure  17.  Circuit  with  R-C  controlled  resonant  frequency. 


It  is  not  the  objective  here  to  exhibit  the  "best"  circuit  for  a  particular  purpose, 
merely  one  that  is  convenient.  Here,  "convenient"  is  taken  to  mean  a  circuit  in  which 
there  is  (at  least  in  principle)  some  independence  in  adjusting  the  b  parameter  and  the 
u-position  of  a  given  pole. 


9.1  R-C  CIRCUIT 

One  such  circuit  is  shown  in  Figure  17.  It  uses  nominal  negative  feedback  which 
becomes  positive  feedback  near  resonance.  The  transmission  of  this  second  order  sys¬ 
tem  from  an  input  V ,  to  outputs  V2,  V3,  and  V4  are 


v2/v,  = 

(s+e)  2/(s:  +2es  +g2) 

(a) 

v3/v,  - 

k  ^  (s+e)/(s2  +2es  +g2) 

0>) 

v4/v,  = 

k1  um2/(s2  +2es  +g2) 

(c) 

(39) 

where  the  various  parameters  in  Eq  (39)  are  given  by 
*  =  A/(l+A) 
u  =  1/RC 

m 

e  =  1/R2C  +  ^/(l+A) 

g2  =  e2  +  ^  K1  (40) 

These  equations  should  be  compared  with  Eqs  (36).  The  obvious  choices  are  to  use 
V3  as  the  output  with  the  odd  pole  configuration,  V  4  with  an  even  pole  configuration. 
However,  each  of  these  outputs  requires  attenuators  to  equalize  an  extra  factor  of  um 
in  Eqs  (40b)  and  (40c). 

With  the  main  poles,  the  preferable  alternative  is  to  use  V  2  for  the  odd  pole  con¬ 
figuration,  V  for  the  even  pole  configuration.  This  leaves  all  of  the  contributions  high 
by  a  factor  of  (s+e).  For  the  main  poles,  a  single  overall  equalization  b/(s  +b)  is  easier 
to  implement  than  individual  equalization  of  the  um.  On  the  other  hand,  the  ampli¬ 
tudes  of  the  corrector  and  compensation  poles  need  adjustment.  Unless  they  have  the 
same  damping  as  the  main  poles,  they  are  best  obtained  from  V  3  (or  V4),  and  equal¬ 
ized  for  combination  with  the  result  of  the  main  pole  synthesis.  The  possibility  of  such 
diversity  in  realization  is  one  of  the  advantages  of  the  pole-residue  approach. 

The  poles  of  the  transmission  given  in  Eqs  (39)  are  at  frequencies  s  m 

s»  =  1/R2C+  V(A+1)±jun  (41) 

To  the  extent  that  the  open  loop  op-amp  gains  A  can  be  made  much  larger  than  the 
ratio  ma/b  at  resonance 


1.  RC  determines  ma. 


2.  R2C  determines  b. 

3.  The  adjustments  of  e  to  b  and  g2-e2  to  m2a2  are  independent  of  one  another. 

Inasmuch  as  all  the  b’s  are  close  to  being  equal,  an  efficient  alignment  procedure  is  to 
fix  R2  and  C,  while  adjusting  R  for  the  individual  frequency.  (This  process  remains 
convenient  for  a  switchcap  equivalent  to  Figure  17.) 

The  amplification  A  is  in  general  frequency  dependent.  Indeed,  it  is  not  uncom¬ 
mon  for  an  op-amp  to  be  in  the  single  pole  roll-off  regime  over  much  of  its  range  of 
use: 


1/(1+A)  =  80  +  s/uu  +  ....higher  order  terms  in  s 

K  =  K0  -s/uu  +  ....  (42) 

where  uu  is  the  frequency  of  projected  unit  gain.  There  are  two  effects. 

1.  An  angle  proportional  to  frequency  is  introduced  into  the  residues. 

2.  There  is  a  reduction  of  the  damping. 

The  former  effect  merely  introduces  a  slight  change  in  the  overall  delay  of  the  filter. 
The  latter  effect  can  be  large  enough  to  cause  oscillation. 

A  stratagem  for  avoiding  this  problem  is  to  use  a  resistance  R  3  in  series  with  C, 
chosen  so  that  R}C  is  l/uu.  This  is  not  a  critical  element.  It  does  not  much  matter 
whether  it  is  physically  in  series  with  C  or  with  the  parallel  combination  of  R  2  and  C. 


9.2  SWITCHCAP  ADAPTATION 

In  integrated  circuits  it  may  be  inconvenient  to  obtain  and  trim  resistors  to  precise 
stable  values.  The  substitute  is  the  switched  capacitor.  Figure  ISa  shows  such  a  capa¬ 
citor  used  to  replace  a  resistor  otherwise  bridged  between  nodes  "A"  and  "B".  In  the 
simplest  case,  the  impedance  looking  out  into  both  nodes  is  low  enough  so  that  the 
settling  time  for  charging  and  discharging  the  switched  capacitor  (C0  in  the  Figure)  is 
short  compared  with  the  time  that  the  switch  is  closed  at  either  node.  The  effective 
conductance  of  the  equivalent  resistor  bridged  between  the  nodes  is  g: 

g  =  C„f  (43) 

where  f  is  the  switching  frequency.  When  a  capacitor  C  intervenes  between  one  termi¬ 
nal  (say,  B)  and  low  impedance  circuitry  (at  say,  D)  then  the  effective  conductance  is 
modified  to 

g  =  C„C  f/(C0+C)  (44) 

The  prime  candidate  for  switchcap  substitutes  are  the  resistors  R  in  Figure  17.  The 
replacement  makes  the  filter  pass  band  programmable  by  adjustment  of  the  switching 
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Figure  18.  Elimination  of  resistors  by  switching  capacitors, 
(a)  Switchcap  circuit,  (b)  Double  switching  circuit. 


rate1  A  second  candidate  for  such  replacement  is  the  damping  control  R  r  Its 
replacement  by  a  switched  capacitor  C,  is  shown  on  the  rhs  of  Figure  18.  Since  there 
are  now  two  capacitors  switched  to  node  B,  it  simplifies  isolation  if  they  are  connected 
to  the  node  alternately  as  shown  in  the  Figure. 

What  are  shown  as  switches  in  the  figure  are,  in  fact,  active  elements.  These  must 
be  designed  so  that  the  settling  time  is  short  compared  with  the  sampling  interval.  In 
addition,  the  closed  loop  response  of  the  overall  system  is  significant  only  "at"  the 
sampling  times,  not  between  them.  As  a  result,  small  delays  in  the  closed  loop  ampli¬ 
fier  do  not  necessarily  show  up  as  additional  delay  in  the  sampled  signal  output. 
Thus,  the  role  of  the  resistors  R  J(  which  is  to  compensate  for  such  small  delays,  may 
be  moot.  In  any  event,  the  precise  value  of  such  resistors  is  not  critical;  if  needed,  they 
can  be  realized  by  any  convenient  means. 


Section  10 

SAMPLED  DATA  FILTERS 


Although  the  topology  of  the  z-plane  filter  has  been  given,  the  analysis  has  been 
carried  out  only  in  the  s-plane.  That  analysis  can  be  extended  to  the  z-plane  topol¬ 
ogy  by  use  of  the  mapping  transformation  given  in  Eq  (9).  This  transformation  maps 
the  z-plane  filter  into  an  infinite  sequence  of  L-Ph  filters  in  the  s-plane.  One  of  these 
is  the  low  pass  filter.  The  remaining  filters  are  "ghosts"  whose  pass  bands  and  singu¬ 
larities  are  2  tt/t  in  (radian)  frequency  apart  along  the  j  u-axis  in  the  s-plane. 

There  are  two  aspects  to  the  resulting  "aliasing"  problem:  (a)  the  cumulative  effect 
of  the  ghosts  within  the  pass-  and  stop-bands  of  the  low  pass  range,  and  (b)  the 
means  by  which  the  transmission  at  the  repeated  "ghost"  pass  bands  is  supressed.  Both 
effects  are  met  in  all  sampled  data  filters,  not  merely  those  which  happen  to  be  L-Ph. 
They  are  treated  here  in  order  to  better  define  the  circumstances  under  which  the 
effects  of  aliasing  are  in  fact  negligible  for  the  L-Ph  filter. 

In  a  band  pass  design,  the  effect  of  one  additional  filter  at  negative  frequencies  is 
neglected.  For  a  low  pass,  the  effect  of  the  next  adjacent  ghost  is  similarly  negligible 
out  to  a  frequency  uc  corresponding  to  one-half  the  sampling  rate: 

we  =  IT. /r  (45) 

If  u0=Ma  is  the  band  edge  of  the  low  pass,  then  the  very  terminology  "low  pass"  pre¬ 
sumes  a  stop  band  wide  compared  with  u0,  so  that 

uc  »  2u0  (46) 

As  a  result  only  the  asymptotic  transmission  of  the  ghosts  affects  the  low  pass  for 

U<UC. 

For  an  odd  pole  number  M,  the  asymptotic  transmission  has  the  form  €b(s-sn)~\ 
where  sn  is  at  the  center  of  the  nth  ghost.  This  asymptotic  form  arises  because  for  odd 
M  the  sum  of  residues  from  all  poles  is  not  automatically  zero.  There  will  be  a  finite 
remainder  «,  whether  by  deliberate  choice  or  by  random  errors  in  adjustment  of  drift 
of  individual  components.  Let  r>( s)  be  the  cumulative  additive  effect  of  the  ghosts. 
Within  the  low  pass  range,  n(jw)  is  mainly  imaginary. 

In  the  pass  band,n(s)  can  be  estimated  by  the  same  methods  used  to  estimate  G(s) 
above.  Near  the  s-plane  origin 

dn/ds  =  €b  /12uc2  (47) 

Since  uc>2u0»  the  effect  in  the  pass-band  cannot  be  greater  than  eb/24Ma,  generally 
less  than  a  fraction  of  a  percent.  Thus  the  ghosts  can  be  neglected  in  the  pass-band. 
For  the  stop-band,  the  effect  of  the  ghosts  is  largest  at  u=uc.  At  ue,  the  effect  can  be 
estimated  by  a  direct  sum  of  the  contributions: 


=  j€b  Z[l/(lnuc-uc)  -l/(2nuc+uc)] 
=  (j«b/uc)  I  l/(n2-l)  n  *  1 
=  .5708  j€b/2Ma 


(48) 


In  this  case,  some  care  must  be  exercised  to  be  sure  that  the  combination  of  €  and  M 
is  appropriate  to  maintaining  the  desired  transmission  floor  in  the  face  of  the  ghost 
contribution. 


A  final  concern  is  the  supression  of  input  frequency  components  which  happen  to 
lie  in  the  ghost  pass-bands.  Ultimately,  a  true  analog  filter  is  required  to  eliminate 
such  components  from  a  non-sampled  data  signal.  However,  it  is  a  filter  whose  band 
edge  can  be  anywhere  between  u0  and  <jc.  No  serious  demands  are  made  on  the  tran¬ 
sition  region  of  this  analog  filter;  it  need  contain  only  a  few  pole  pairs. 


The  switchcap  technique  for  realizing  2nd  order  operations  is  a  digital-analog 
hybrid  in  which  the  samples  are  held  between  sampling  instants.  The  result  is  to  mul¬ 
tiply  the  periodic  s-plane  spectrum  by  a  sinx/x  whose  zeros  are  exactly  at  the  centers 
of  the  ghost  pass-bands  along  the  ju-axis.10 


At  the  edges  of  the  ghost  bands,  the  transmission  is  bounded  by  that  of  the  closest 
edge  of  the  nearest  ghost  to  Tgho<1 


1  ghott 


=  TT  Up  /2(J 


(49) 


uc  is  arbi- 


Tgh  can  be  made  arbitrarily  small  by  raising  the  sampling  rate  so  that 
trarily  high.  However,  the  sampling  rate  may  be  fixed  by  practical  considerations;  in 
such  cases  the  additional  analog  filter  at  the  input  is  inavoidable. 
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10  With  a  periodic  spectrum  this  multiplication  by  sinx/x  does  not  produce  divergence 
as  s-±°°;  the  periodic  filter  response  decreases  exponentially  in  this  limit  more 
strongly  than  the  sinho/o  increases. 


Section  11 


DISCUSSION  OF  RESULTS 


The  question  most  often  asked  about  the  Linear  Phase  Filter  is,  "how  does  it  com¬ 
pare  with  a  standard  Butterworth  or  Tchebychef  design  plus  all-pass  phase  equaliza¬ 
tion?" 


The  comparisons  made  in  the  text  between  the  simple  L-Ph  filter  and 
the  Tchebychef,  and  also  between  the  feedback  filter  and  the  Butterworth, 
should  suffice  for  the  phase  equalization  comparison.  The  comparison  is 
most  favorable  when  the  tolerable  amplitude  and  phase  errors  over  most 
of  the  pass  band  are  small,  so  that  Tchebychef,  Butterworth,  and  L-Ph 
design  procedures  yield  nearly  the  same  results  for  the  first  SO  dB  or  so  of 
stop  band  attenuation.  For  such  cases,  the  L-Ph  filter  needs  no  phase 
equalization  in  the  pass  band,  yet  has  as  few,  or  fewer,  second  order  ele¬ 
ments  than  the  unequalized  competition. 

Another  questions  is,  "what  are  the  disadvantages  of  these  filters?"  The  answer 
depends  on  the  application  and  the  components  to  be  used. 


One  potential  disadvantage  of  the  L-Ph  designs  is  that  they  all  achieve 
part  of  their  stop  band  attenuation  by  balance  between  YA  and  YB. 
These  two  transmissions  must  be  made  of  components  that  drift  alike  with 
temperature  and  time.  In  particular,  they  should  be  mounted  on  the  same 
chip  symmetrically  with  respect  to  sources  of  heat. 

Note  that  the  figures  in  the  text  are  all  normalized  to  unity  transmis¬ 
sion  at  d.c.  The  absolute  transmissions  of  the  simple  filters  are  low,  a 
consequence  of  b/a^l  as  shown  in  Figure  8.  As  a  result,  these  filters  are 
more  sensitive  to  drift  than  the  loop  topologies  for  which  b/a  can  be 
taken  as  close  to  zero  as  the  stability  of  components  permits. 

Another  potential  disadvantage  is  that  most  L-Ph  designs  have  floors 
(or  at  least  landings  on  a  descending  ramp)  to  the  stop  band  transmission. 
It  is  important  to  recognize  that  such  floors  also  arise  from  what  is  prac¬ 
tical  in  the  isolation  of  input  from  output  in  a  single  filter  chip  or  box, 
regardless  of  theoretical  transmission.  Such  isolation  floors  range  from  as 
little  as  40  dB  to  as  much  as  70  dB,  depending  on  the  frequency  range 
and  the  technology  available.  The  fact  that  individual  L-Ph  filters  have 
theoretical  floors  or  landings  in  the  same  range  cannot  be  counted  as  a 
severe  additional  disadvantage.  The  cure  is  in  both  cases  the  same:  if 
more  isolation  is  necessary,  use  two  stages  in  cascade. 


It  is  sometimes  asked,  "how  well  do  these  filters  conform  to  a  suitable  set  of  theo¬ 
retical  limitations?" 

Such  questions  can  be  answered  only  in  the  context  of  additional  intellec¬ 
tual  assumptions  or  esthetic  standards.  A  framework  for  such  discussions 
was  provided  in  the  subsection  on  multiple  pole  correctors  in  the  form  of 
a  transition  zone  within  which  the  transmission  is  only  loosely  controlled. 

However,  a  detailed  discussion  of  the  topic  is  beyond  the  scope  of  the 
present  study. 

A  final  question  is  "why  aren’t  these  L-Ph  designs  more  widely  used?" 

The  answer  is  essentially  beaurocratic.  Filter  specifications  are  written 
with  other  filters  in  mind;  enough  features  of  such  other  filters  are  speci¬ 
fied  to  make  the  use  of  any  other  filter,  including  the  L-Ph  filter,  awk¬ 
ward. 

The  problem  of  specification  is  two-fold,  amplitude  and  phase.  The  amplitude  specifica¬ 
tion  problem  is  merely  one  of  making  engineers  aware  of  the  fact  that  for  a  small  adjustment 
of  their  transmission  template,  they  can  get  L-Ph  "for  free". 

Most  of  the  examples  given  in  this  report  easily  fit  the  amplitude  template  shown  in  Fig¬ 
ure  1.  This  template  differs  from  the  template  that  might  accompany  a  Butterworth  or 
Tchebychef  design  in 

a)  Specifying  a  1  dB  (5‘  phase  distortion)  point  instead  of  a  3  dB  point. 

b)  Allotting  more  space  between  the  in-band  limits  and  the  1  dB  point  than  is  typical  of 
Butterworth  or  Tchebychef  designs. 

c)  Allowing  an  isolated  deep  stop  band  zero  (carrier  supression)  up  to  the  edge  of  the 
transition  zone. 

d)  Not  specifying  stop  band  attenuation  deeper  than  fixed  floor  from  a  single  stage  fil¬ 
ter. 

From  several  of  the  specific  designs  that  appear  in  the  text,  it  is  evident  that  the  template 
is  conservative.  The  floor  on  the  progression  of  band-edge  carrier-supression  zeros  is  based 
on  nominal  damping  for  a  nominal  out-of-band  zero.  If  further  in-band  phase  distortion  is 
tolerable,  or  if  higher  Q’s  are  practical  for  the  zero,  then  this  floor  is  lowered. 

The  phase  template  for  these  filters  is  the  same  as  the  amplitude  template,  except  that 
the  1  dB  point  is  replaced  by  the  5*  or  0.1  radian  point.  There  is  also  no  specification  of 
the  phase  beyond  this  point.  The  specification  of  phase  is  made  difficult  by  the  common  use 
of  "group  delay"  defined  as  the  frequency  derivative  of  phase.  Regrettably,  "group  delay"  so 
defined  has  no  physical  significance  when  the  phase  distortion  is  low.  11  Yet  group  delay  (by 
derivative)  distortion  is  almost  always  specified. 


Group  delay  acquires  physical  meaning  when  the  Fourier  transform  for  the  impulse 
response  can  be  found  by  the  method  of  stationary  phase  integration.  This  requires 
that  the  phase  distortion  be  high,  amounting  to  at  least  several  radians  across  the 
frequency  bandwidth  of  the  "group"  being  considered. 


ii 


These  Alters  behave  in  the  pass-band  like  a  piece  of  cable  with  a  Axed  delay  and  small 
impedance  mismatches  at  the  two  ends.  In  the  simple  topology,  the  effective  mismatch  is 
inherent,  producing  the  periodic  errors  discussed  in  Section  4.  In  the  loop  topologies,  the 
transmission  line  is  theoretically  terminated  at  one  (single  loop)  or  both  (two  loop)  end(s). 
In  the  loop  topologies,  mis-settings  of  gains  G0  will  cause  the  mismatches  to  return.  As  a 
result  of  the  mismatches,  there  are  internal  reAections  which  are  small,  broadband,  and  cir¬ 
culate  with  decreasing  amplitude  at  each  bounce.  The  delay  of  the  primary  echo  is  simply 
triple  travel;  higher  order  echos  have  correspondingly  longer  multiple  travel  delays. 

How  big  an  echo  can  one  tolerate?  That  depends  on  the  application.  In  an  audio  chan¬ 
nel  a  broad  band  30  ms  echo  down  6  dB  is  scarcely  distinguishable  from  the  environment  by 
the  ear,  but  is  disastrous  to  digital  data  transmission.  On  the  other  hand  a  1  s  echo  down 
30  dB  can  be  disastrous  to  listening  (or  speaking)  and  scarcely  noticable  in  its  effect  on 
digital  data  transmission.  To  multiply  the  1  s  echo  by  its  amplitude  and  announce  a  30  ms 
delay  distortion  is  not  only  absurd:  it  leads  to  the  wrong  engineering  conclusions.  Yet  this  is 
exactly  what  a  "group  delay  distortion"  calculation  does. 

A  solution  to  the  phase  distortion  speciAcation  problem  is  to  specify  the  phase  distortion 
in  two  parts 

a)  Recognizable  broad  band  echos,  and 

b)  Other  sources  of  delay  distortion. 

Examination  of  the  Figures  in  this  report  reveals  that  the  pass  band  amplitude  and  phase 
ripple  are  strongly  correlated.  This  is  because  both  consist  almost  entirely  of  broadband 
echos.  Thus,  the  phase  distortion  can  be  speciAed  by  echo  amplitude  or  by  radian  error;  but 
the  amplitude  of  the  identiAable  echo  can  be  frequency  dependent.  The  "other  sources"  of 
delay  distortion  are  not  signiAcant  to  within  a  units  of  the  in-band  zone  edge. 

It  is  now  more  than  thirty  years  since  the  writer  built  his  Arst  L-Ph  band-pass  and  low- 
pass  Alters.  The  statements  made  in  his  1964  paper  are  still  valid:  if  there  is  a  use  for  these 
designs,  they  will  be  used;  if  not,  they  won’t  be.  The  only  valid  concern  is  that  these  Alters 
be  a  known  alternative,  so  that  the  design  engineer  in  fact  has  a  choice. 
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